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Abstract. 

We show that terrestrial planets could survive in variously restricted 
regions of the habitable zones of 47 Ursae Majoris, e Eridani, and p 
Coronae Borealis, but nowhere in the habitable zones of Gliese 876 and v 
Andromedae. The first three systems between them are representative of 
a large proportion of the 90 or so extrasolar planetary systems discovered 
by mid-2002, and thus there are many known systems worth searching 
for terrestrial planets in habitable zones. We reach our conclusions by 
launching putative Earth-mass planets in various orbits and following 
their fate with a mixed-variable symplectic integrator. 


1. Introduction 

It will be some years before we know whether any of the known extrasolar 
planetary systems have planets with masses the order of that of the Earth. We 
have therefore used a computer model to investigate whether such planets could 
be present, and in particular whether they could remain confined to the habitable 
zones (HZ) of a representative sample of the known extrasolar planetary systems. 
If so then it is possible that life is present on any such planets. 

The HZ is that range of distances from a star where water at the surface of 
an Earth-like planet would be in the liquid phase. We have used the intermediate 
boundaries for the HZ originating with Kasting, Whitmire, Sz Reynolds (1993). 
Because of simplifications in the climate model, these distances are conservative 
in that the HZ is likely to be wider. For zero-age main-sequence stars (ZAMS 
stars) the HZ lies closer to the star the later its spectral type. As the star ages 
the boundaries move outwards. 

To test for confinement to the HZ, we launch putative terrestrial planets, 
labelled EM, into various orbits in the HZ and use a mixed-variable symplectic 
integrator (Chambers 1999) to calculate the evolution of the orbit. The integra¬ 
tion is halted when a terrestrial planet comes within three Hill radii (3i?n) of the 
giant. This is the distance at which a symplectic integrator becomes inaccurate, 
and is also the distance at which severe orbital perturbation of the terrestrial 
planet occurs. Confinement requires that the integration lasts for the full pre-set 
time, and that the semimajor axis of the terrestrial planet remains within the 
HZ. The pre-set time depends on the shortest orbital period in the system. The 
integrator requires a time step shorter than about 5% of this shortest period, 
and the pre-set time must keep the CPU time per integration less than 100 or 
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so hours. For systems like p Coronae Borealis the orbital period of the giant 
planet is so short that we use 100 Ma (100 million years) as the standard, and 
also for Gliese 876 and v Andromedae. For 47 Ursae Majoris and e Eridani we 
use 1000 Ma as the standard. 


2. The exoplanetary systems studied 

We selected five systems that between them represent a large proportion of 
the 90 or so extrasolar planetary systems known by mid-2002. Figure 1 shows 
the important characteristics of these systems. All five systems have been seen 
only by Doppler spectroscopy, which yields m sin(?'o) rather than m. For an 
integration we put in an actual mass m for the giant planet. This is tantamount 
to setting iq to some value. We used = 90°, which gives the minimum value 
for m, and other values that varied from one system to another. 

All planets are launched with zero mean anomalies but with various lon¬ 
gitudes of the periastra. The outcome can be sensitive to the differences A w 
between these longitudes, in that some differences result in early close encoun¬ 
ters (EM comes within 3i?n) of the giant planet, whereas other differences result 
in no close encounters within the pre-set integration time. 


3. Results 

For Gliese 876 and v And inspection of Figure 1 shows that the whole of HZ(0) 
and HZ (now) are traversed by (3 Ru + A r), and so there should be no chance 
of avoiding an early close encounter between EM and a giant planet. This 
expectation is born out by the integrations, with close encounters no later than 
about 10 Ma. The second giant now known around Gliese 876 makes things 
worse. 

For p CrB at minimum giant planet mass confined orbits for EM are found 
right across HZ(0), even at the 1:6 mean-motion resonance near the inner bound¬ 
ary of HZ(0). HZ(now) is further out and so here too confinement is ubiquitous. 
We obtain similar confinement at eight times the minimum giant mass. 

The age of p CrB is about 6000 Ma, so confinement for the standard pre¬ 
set integration time of 100 Ma could be a poor indicator of confinement for 
6000 Ma. We therefore extended a few integrations for pre-set times 2-3 times 
longer than the standard. Full-term confinement was obtained in every case. 
Moreover, there is no sign of any long-term secular increase in the eccentricity 
of EM. Therefore, we believe that confinement for 6000 Ma is common, even 
ubiquitous, in the habitable zones of p CrB. 

e Eri is a K2V star and therefore evolves more slowly than the Sun, and it is 
young, only 500-1000 Ma. Therefore HZ(0) ~ HZ(now). We explored the inner 
part of the habitable zone and the region immediately interior to it, including all 
the mean motion resonances, which are closely spaced. At such resonances there 
is a strong tendency for A®emg = 180° to give confinement, and for Azuemg 
= 0° to fail to do so, a well-known result. 

At minimum giant planet mass, EM orbits that remain confined for the 
1000 Ma pre-set time, are restricted to initial (t = 0) semimajor axes of EM, 



Terrestrial Planets in Known Extrasolar Planetary Systems 


3 


Rho CrB >1.1 Mj 

0.72 

0.94| 

1.42 

] 1.85 

_ 1 _ 


i 

i 

i 

_i_ i _ 

_ 1 _ 

i 

i 

i 


0 0.5 distance/AU 1.0 1.5 2.0 


0.47 

l 

0.91 s Eridani > 0.86 M [ 

i - m - 


1 1 

0 

i 

2 3 4 

distance/AU 

5 6 

47 UMa 

0.76 

i.ooj 

1 

1.43 11.90 

i - 




1 

_ 1 _ 

T , > 2.54 M, 

> 0.76 M s 


0 1.0 distance/AU 2.0 3.0 4.0 


Gliese 876 

0.112 

0.221 

> 2.1 M ' 

1 _ 

-#— 1- 

r 

i 

_i_ i _ 

_*_ T 

1 

_ 1 _ i _ 1 _ 


0 0.1 distance/AU 0.2 0.3 0.4 


Ups And 1-25 

• - 0 - 

1.80 

i 

2.49 

-•- 

I 3.50 

_L 

>0.71 Mj >2.11 Mj 

1 

> 4.61 Mj 

1 


0 1.0 distance/AU 2.0 3.0 4.0 


Figure 1. Characteristics of the five exoplanetary systems studied. 
The ZAMS habitable zone HZ(0) is shown shaded, and the boundaries 
of HZ (now) by vertical dashed lines. Each giant planet is shown by 
a black disk labeled with the value of m sin(io), where m is the mass 
of the giant planet and iq is the inclination of the planet’s orbit with 
respect to the plane of the sky. At each giant planet the solid line shows 
the total excursion 2Ar of the giant planet due to its eccentricity. The 
dashed line extends to (3 Rh + Ar) each side of the giant planet when 
it has its minimum mass (Ru is proportional to m 1//3 ). 
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«em(0) < 0.59 AU. Within this distance there is an intricate pattern of outcomes, 
depending on whether Azuemg = 180° or 0°, and whether oem(0) is at a mean- 
motion resonance. Only within 0.44 AU is there general confinement. We also 
explored 1.39 times the minimum mass. This is because e Eri has a circumstellar 
dust ring that is seen to be elliptical. Assuming it to be circular we get an 
estimate of ?'o = 46°, that corresponds to the 1.39 multiplier. At this greater 
mass there is general confinement within 0.44 AU, but no confinement at all at 
&em(0) > 0.46 AU. 

47 UMa has an outer giant with a present eccentricity e < 0.2. In Figure 
1 the value is 0.1. This value was also used in the integrations at t = 0, and 
also a value of 10 -5 . As well as minimum mass, the giant planets were given 
1.5 times the minimum, in accord with estimates of zq for this system (Gonzalez 
1998). At 1.5 minimum mass, with a few exceptions, the two giant planets were 
launched with AtUGiG 2 = 0° and 120°, and Axuemgi = 120° and 180°. This 
taught us that AotgiG 2 = 120° is less likely to lead to confinement. So, to err 
on the side of caution in our claims for confinement, at minimum mass we ran 
configurations other than 120° only if there was no confinement at 120°. 

For both masses there is no confinement beyond 1.25 AU. At minimum 
mass, at most values of oem(0) < 1.25 AU we obtain confinement at both 
eccentricities and at all values of the A w. At 1.5 times the minimum, with 
e(0) = 10 -5 the outcome is similar at uem(0) < 1.20 AU, but at e(0) = 0.1 
confinement is obtained only at a few values of uem(0) < 1.15 AU and for some 
Aw. 

47 UMa is about 7000 Ma old (Gonzalez 1998). The pre-set integration 
time of 1000 Ma is a modest fraction of this, though we see close encounters 
within 500 Ma in all but a small proportion of cases, and in these cases there 
are clear upward trends in eccentricity of EM. Therefore, confinement for 1000 
Ma is very likely to mean confinement for 7000 Ma. 

Further details of our work on all but the e Eri system are at Jones, Sleep, 
& Chambers (2001) and Jones & Sleep (2002). 


4. Conclusions 

Confinement of a terrestrial planet for the age of the star is likely throughout 
the ZAMS and present-day habitable zones of p CrB. In e Eri (where the habit¬ 
able zone has barely moved outwards) only the innermost part of the HZ offers 
confinement at minimum giant planet mass, and then only for certain launch 
configurations of the terrestrial body. At 1.39 times this minimum mass, con¬ 
finement is restricted to the inner boundary of the HZ (and inwards), again for 
certain configurations only. In 47 UMa confinement is found in the inner part of 
the present HZ, particularly if either the present eccentricity is much less than 
0.1 or if the giant planets’ masses are close to the minimum values, or if both 
conditions hold. Gliese 876 and v And are devoid of confined orbits in their 
habitable zones. 

Our conclusions on 47 UMa are in accord with work done by others, though 
this other work is significantly less extensive (Laughlin, Chambers, & Fischer 
2002; Noble, Musielak, & Cuntz 2002). 
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